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FOREWORD 


The  study  presented  in  this  report  covers  a  phase  of  a 
theoretical  investigation  of  curved  jets  supported  by  the  U.  S.  Army 
transportation  Corps  during  the  year  1960-lii61  through  the  Air 
Programs  Office  of  the  Office  of  Naval  Research,  under  the  terms 
of  Contract  Nonr-2747(00). 
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SUMMARY 


A  numerical  solution  of  the  two-dimensional  curved  turbulent 
and  incompressible  jet  flow,  using  the  IBM  7090  computer,  is  presented. 
This  solution  is  a  straight-forward  extension  of  the  classical  straight 
jet  solution.  Its  only  limitation  is  that  it  assumes  similar  velocity 
profiles. 

The  solution  is  obtained  by  reducing  the  Navier  Stokes  partial 
differential  equations  of  the  curved  jet  flow  to  a  third-order  total  different¬ 
ial  equation  of  the  Hartree-Skan  type.  This  equation  is  integrated  numeric¬ 
ally  on  the  IBM  7090  computer  using  a  Runge-Kutta  subroutine.  The 
boundary  conditions  are  carefully  established  and  discussed.  The  complete 
program  for  the  numerical  solution  is  included  in  the  report.  A  duplicate 
IBM  deck  can  be  obtained  from  Aerophysics  Company  upon  request. 

Typical  numerical  results,  in  the  form  of  velocity  and  pressure  distri¬ 
bution  across  the  jet  are  presented  and  discussed. 
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maximum  jet  velocity 

u 

jet  velocity  component  in  x-direction 

V 
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*  1 

curvilinear  system  of  coordinates,  as  defined  in  the  text  and 
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on  figure  1 

a 
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INTRODUCTION 


In  spite  of  the  remarkable  success  of  inviscid  theories  to 
predict  the  performance  of  hovering  ground  effect  machines,  it  cannot 
be  concluded  that  viscous  phenomena  do  not  sometimes  play  an  important 
role  in  the  understanding  of  the  flow  mechanics  of  hovering  annular  jets 
in  ground  effect.  In  particular,  the  lack  of  an  analytical  grasp  of  viscous 
effects  has  prevented  until  now  the  establishment  of  a  satisfactory 
mathematical  theory  of  the  stability  of  ground  effect  machines. 

The  major  viscous  effect  associated  with  annular  jet  flow  concerns 
entrainment  of  air  from  the  central  cavity,  the  so-called  "air  cushion", 
by  the  peripheral  jet.  When  the  base  of  the  machine  is  parallel  to  the 
ground,  this  entrainment  results  in  a  quasi -rotational  flow  under  the 
machine,  which  can  be  ignored  for  most  practical  purposes,  especially 
if  the  base  is  fairly  close  to  the  ground.  However,  when  the  base  of  the 
machine  is  pitched  with  respect  to  the  ground,  a  strong  cross-flow 
establishes  itself,  with  an  unsymmetrical  pattern,  and  strong  destabil¬ 
izing  pitching  and  rolling  moments  may  appear.  The  seriousness  of 
this  phenomenon  depends  obviously  upon  the  amount  of  air  which  is  being 
entrained  by  the  peripheral  annular  jet. 

In  a  very  elegant  first  approximation,  Chaplin  assumed 
(Reference  1)  that  the  amount  of  flow  entrained  by  a  curved  jet  was  the 
same  as  that  of  a  straight  jet,  for  which  both  theoretical  and  experiment¬ 
al  data  are  available  (References  2  and  3).  The  purpose  of  the  present 
investigation  is  to  approach  the  problem  analytically,  by  extendit^  the 
classical  straight  jet  flow  theory  to  the  curved  jet. 

A  formulation  of  the  problem  was  made  earlier  by  Boehler  in 
Reference  4  for  a  laminar  jet  and  by  Chang  in  Reference  5  for  a  turbulent 
jet.  Here,  the  turbulent  case  is  considered,  since  it  is  well  known  that 


all  jet  flows  obtained  experimentally  are  mostly  turbulent.  Formally, 
the  equations  of  the  present  problem  are  the  same  as  those  of  Reference 
5.  However,  an  important  difference  with  Reference  5  is  that  the  systems 
of  axes  are  defined  differently.  In  all  cases,  the  basic  assumption  is 
made  of  the  similarity  of  the  velocity  profiles.  This  assumption  is 
restrictive,  and  may  turn  out  in  the  long  run  not  to  be  justified. 

However,  it  is  a  reasonable  one  for  a  first  approach  to  the  problem, 
since  it  makes  possible  a  closed-form  solution  of  the  problem  and 
allows  a  simple  direct  comparison  with  the  straight  jet  case.  Restrictive 
as  it  may  be,  the  assumption  of  similarity  was  also  used  in  the  analytical 
solution  of  the  boundary  layer  along  a  curved  wall  (Reference  6). 

The  solution  is  obtained  by  reducing  the  Navier-Stokes  partial 
differential  equations  of  the  curved  jet  flow  to  a  third-order  total 
differential  equation  of  the  Hartree-Skan  type.  This  equation  is  inte¬ 
grated  numerically  on  the  IBM  70:^0  computer  using  a  Runge-Kutta 
subroutine.  The  boundary  conditions  are  carefully  established  and 
discussed.  The  complete  program  for  the  numerical  solution  is 
included  in  the  report.  Typical  numerical  results,  in  the  form  of 
velocity  and  pressure  distribution  across  the  jet  are  presented  and 
discussed. 

The  solution  of  the  basic  curved  viscous  jet  is  important  not 
only  for  an  understanding  of  the  flow  mechanics  of  ground  effect  machines, 
but  for  many  other  problems  as  well,  such  as  those  involving  cavity  flows 
(base  pressure  of  missile  boosters  for  example)  or  the  understanding 
of  the  Coanda  effect. 

The  two-dimensional  solution  can  be  extended  to  three  dimensions, 
if  desired,  by  means  of  a  Mangier  transformation  (Reference  4). 
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I  EQUATION  OF  THE  PROBLEM 


Historically,  it  was  found  (Reference  2,  page  143  for  example) 
that  the  equations  governing  the  flow  of  a  two-dimensional  jet  emerging 
from  a  narrow  slit  into  a  fluid  at  rest  are  formally  the  same  as  the 
boundary  layer  equations  along  a  straight  wall  (with  no  pressure 
gradient);  actually  the  two-dimensional  jet  case  is  one  of  the  few 
cases  where  the  Navier-Stokes  equations  can  be  integrated  exactly. 

Similarly,  the  equations  of  motion  of  a  two-dimensional  curved 
viscous  jet  are  formally  the  same  as  those  of  the  boundary  layer  along 
a  curved  wall.  The  general  Navier-Stokes  equations  for  the  flow  along  a 
curved  wall  and  their  simplification  by  use  of  the  standard  boundary- 
layer  approximations  are  discussed  in  detail  in  Reference  6.  These 
equations  are  directly  usable  for  the  two-dimensional  curved  jet,  as 
pointed  out  in  Reference  4.  One  difference,  however,  lies  in  the  defini¬ 
tion  of  the  system  of  axes  and  coordinates.  This  point  bears  a  detailed 
discussion  which  will  be  made  here. 


For  the  boundary-layer  flow  along  a  curved  wall,  the  choice  of 
the  coordinate  system  is  obvious  (figure  la).  It  is  a  curvilinear  system 
of  coordinates  in  which  the  x  axis  is  in  the  direction  of  the  wall  and  the  y 
axis  is  normal  to  it.  The  stagnation  point  is  taken  as  origin  of  coordinates. 
The  curvilinear  net  therefore  consists  of  curves  parallel  to  the  wall  and 
of  straight  lines  perpendicular  to  them.  The  corresponding  velocity 
components  are  called  u  and  v  and  the  curvature  at  a  point  x  is  K(x), 
which  is  a  continuous  function  of  x  and  is  positive  for  convex  and 
negative  for  concave  curvature.  For  the  laminar  case  and  for  moderate 
wall  curvatures,  it  is  shown  in  Reference  6  that  the  boundary  layer 
equations  reduce  to: 
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Similarly,  for  the  turbulent  case,  the  equations  of  motion 
reduce  to: 


^  u  Du  ^  V  "Ju 
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€  is  the  kinematic  viscosity  which  is  assumed  to  be  determined 
from  Prandtl's  mixing-length  theory  (Reference  2,  Chapter  XIX).  In 
general: 


€  *  X|  b  -  UjniP 

here:  (IH)  f -X,  b  U 


where  U  denotes  the  maximum  velocity,  b  the  width  of  the  mixing  zone 
and  X,is  an  empirical  constant.  C  is  assumed  to  be  independent  of  y. 


Denoting  the  center-line  velocity  and  the  width  of  the  jet  at  a 
fixed  characteristic  distance  s  from  the  orifice  by  U„  and  b„  respective- 

9  9 

ly,  we  may  write: 


U  *  U 


s 


b  ■  b„  X 

®  -s 


Consequently: 

^  -  ^s 


with 


Further,  the  assumption  of  similarity  results  in  the  turbulent  case  in: 

(IV)  7?  -O'JL 

X 


where 


^  is  a  free  constant. 
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We  may  now  come  back  to  the  problem  at  hand,  i.  e.  the  curved 
jet  (fig.  lb).  The  main  difficulty  in  trying  to  apply  the  above  boundary 
layer  equations  lies  in  the  definition  of  a  system  of  coordinates.  One  is 
first  led  to  separate  the  jet  into  two  regions:  the  "Inside"  jet  and  the 
"outside"  jet.  The  line  which  separates  the  two  r^ions  is  then  defined 
as  the  "axis"  of  the  jet.  At  the  outset,  the  shape  of  the  axis  is  not  known. 
It  will  be  shown  later  on  that,  out  of  many  possible  such  axes,  only  one 
must  be 'chosen,  for  no  better  reason  than  because  that  particular  choice 
leads  to  a  possible  reduction  of  the  Navier-Stokes  partial  differential 
equations  to  a  single  total  differential  equation.  Thus  a  two-dimensional 
curved  jet  flow  is  the  sum  of  two  flows  respectively  analogous  to  the 
boundary-layer  flow  along  a  concave  wall  of  curvature  -  R  and  that  along 
a  convex  wall  of  curvature  -t-  R.  For  each  region,  then,  using  the  proper 
sign  of  R,  the  above  equations  I  and  II  hold. 

One  would  be  tempted  to  say  that  the  maximum  velocity  of  the 
flow  U  defined  before  in  connection  with  the  definition  of  the  kinematic 
viscosity  (equation  EQ)  is  the  center-line  velocity,  i.  e.  the  velocity  along 
the  "axis"  of  the  jet,  as  is  the  case  for  the  straight  jet.  We  shall  see 
later  that  this  is  the  case.  Indeed,  one  could  arbitrarily  define  the 
axis  as  the  line  of  maximum  velocity;  however,  this  condition  would  be 
insufficient  to  determine  the  equation  of  the  axis  and  would  lead  to 
analytical  difficulties  in  the  solution  of  the  problem. 

By  analogy  with  equation  (IV)  of  the  boundary  layer  case,  the 
assumption  of  similarity  of  the  velocity  profiles  is  used  in  the  solution 
of  the  curved  jet  case.  The  hypotheses  of  Prandtl's  mixing-length 
theory  are  also  assumed  to  hold,  possibly  with  different  values  of  the 
arbitrary  constants,  so  that  the  maximum  velocity  U  still  varies  like 
the  inverse  of  the  square  root  of  x. 

In  the  algebraic  calculations  involved  in  passing  from  the  Navier- 
Stokes  partial  differential  equations  to  the  total  differential  equation,  it  is 
sufficient  to  consider  only  one-half  of  the  problem,  i.  e.  either  the 
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inside  or  the  outside  jet,  since  the  other  half  can  be  obtained  by  changing 
the  sign  of  the  radius  of  curvature  in  the  equations. 

The  problem  of  the  curved  viscous  jet  can  therefore  be  stated  as 

follows: 


■Ju 
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To  solve  these  equations,  one  can,  as  for  the  straight  jet,  define 


a  stream  function, 

,  1 

y.  U3  s’* 


which  is  assumed  to  be  of  the  form: 


1 

F(7|) 


(6) 


F(7|)  is  the  unknown  function,  the  determination  of  which  will  give  the 
solution  of  the  problem. 


Expressing  the  velocity  components  in  terms  of  F  and  F'  and 
substituting  into  the  momentum  equations,  one  can  replace  the  partial 
differential  equations  by  one  total  differential  equation  involving  the  function 
F.  It  must  be  noted  that,  in  the  following  derivations,  some  liberties  have 
been  taken  with  respect  to  the  notations  for  partial  differentiation.  These 
equations  have  been  checked  very  carefully,  however,  and  are  believed 
to  be  correct.  The  calculations  are  as  follows: 

■  Ug(_s_)i  F(7|  ) 
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^y  *5  y 
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(7) 


T  X 


U 
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(8) 


from  equations  (T^and  (8),  the  left-hand  member  of  equation  (1)  can 
easily  be  calculated. 

Next  one  must  calculate ,  using  equation  (2)  and  substitute 
the  value  thus  found  into  the  right.-hand  member  of  equation  (1).  In  keeping 
with  the  general  assumption  of  similar  velocity  profiles,  one  will  assume  that 
y  only  appears  in  the  expression  of  p  through  the  term7|  .  Hence: 


29  21 


Dx 


Further: 
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T>y  n 
From  equation  (2) : 


(  10) 
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From  equation  (7): 

1 

u-Us(4)^  F'(7^)  (11) 

We  stated  earlier  that  y  could  not  appear  formally  in  the 
expression  of  p;  similarly,  it  cannot  appear  in  the  e}q;)re8sio  n 
of -^8^.  Therefore,  equation  (11)  is  only  meaningful  if  the  product  Ky  is 
a  function  of  7^  alone.  Rather  than  attempting  to  solve  the  problem  for 
any  function  of  7^  ,  one  can  choose  a  particular  one,  the  simplest  one,  for 
example  a  linear  function  of  .  One  is  therefore  led  to  put: 


Ky 


4^? 


(  12) 


where  A  is  a  constant. 
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or: 


2  Ky-  A  ri 


(12b) 


By  substitution,  equation  (11)  simplifies  to; 

Tf*  (13) 

Note  that,  now,  equation  (12)  defines  the  center  line  of  the  jet. 

By  substitution,  one  finds  the  following  differential  equation  for  F: 
(A  +  1).  F’2  +  F  F"  +  I  F’”  »  0  ^  (  14  ) 

after  letting: 

(  15) 

which  is  legitimate  since  X,  in  the  expression  of  €  ^  is  a  free  constant. 

Naturally,  from  equation  (12),  if  A»  0,  the  curvature  K  of  the  jet 
is  also  zero.  A  «  0  corresponds  to  the  straight  jet  case.  Equation  (14) 
reduces  then  to: 

f2  +  F  F"  +  ^  F"’  -  0  (  16  ) 

Equation  (16)  is  the  differential  equation  of  the  straight  jet  solution  (Ref¬ 
erence  2,  page  4ii8),  as  it  should  be. 

The  solution  of  the  curved  jet  problem  amounts  to  finding  a  solu¬ 
tion  of  equation  (14),  which  also  satisfies  the  proper  boundary  conditions. 
This  solution  must  be  referred  to  a  system  of  axes  which  still  remains  to 
be  determined.  What  one  knows  at  this  point  is  that  the  axes  are  defined 
by  equation  (12  b):  2  Ky  «  ATJ  .  Before  proceedir^  with  the  solution, 
we  shall  determine  the  center  line  of  the  jet  and  these  axes  analytically. 
Obviously,  the  center  line  will  consist  of  a  family  of  curves  depending  upon 
one  arbitrary  parameter  A. 
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n  SHAPE  OF  THE  JET 


The  shape  of  the  jet  is  defined  by  the  equation 

Ky.^7) 

where:  x  and  y  are  curvilinear  coordinates  along  the  curve,  as  shown 
on  the  figure 

K  is  the  curvature  of  the  center  line  of  the  jet 
l\m  y/x,  where  Cf  is  a  constant 
A  is  a  parameter 


When  one  looks  for  the  equation  of  the  shape  of  the  center  line  of  the  jet, 
which  we  have  called  x-axis,  one  must  find  this  shape  with  respect  to 
a  fixed  system  of  rectangular  coordinates  x^  and  y^*  The  unknown 
equation  of  the  curve  is  y^  «  y^  (x)  and  is  defined  by  the  following  intrin¬ 
sic  equation  (since  x  becomes  s^,  length  of  an  arc  on  the  curve): 

Ky  ■  A/2  O'  (y/sj) 
or  K  •  A  O'  /2  (l/sj) 

Defining  a  new  constant  C  «  A  (S / 2,  one  has: 


K  *  C/sj 
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Therefore,  one  wishes  to  find  the  curves  for  which  the  curvature 
at  a  given  point  is  inversely  proportional  to  the  length  of  the  arc  between 
the  origin  and  that  point,  the  length  of*  the  arc  being  defined  by  the  differ¬ 
ential  relationship; 

dsj^  ■  dxj^  +  dyj^ 


For  simplification,  we  shall  drop  the  suffix  "1”  in  all  that  follows.  Let 
us  remember,  though,  that  the  curve  to  be  found  will  be  referred  to  the 
fixed  system  of  coordinates  x^Oy^. 


The  radius  of  curvature  R 
o  3/2 


1/K,  is  given  by 


The  intrinsic  equation  is: 
K*  C/s 


(  17) 


(  18) 


Therefore, 

R«  (1/C)s 

Differentiating  with  respect  to  x,  there  comes: 

dR  .  1  ds 
cGT  *C1SBr 

but  ^  ^  VTTT^ 

dx 

and  dR  _  3/2  [l  +  y’2]  ^  x  2y'y'’^  -  y’”  [  1  +  y’^1 

^  . - - - - 

r  jl/2  o  r  ol3/2 

or:  y  ^  3  y*2|  x  yY’2  -  y’»'  \  1  y’2j _ 

dx  ^tf2 


( lij) 

(20) 

(21) 

(22) 

(23) 
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Substituting  (21^  and  (2|  into  equation  (4),  there  comes: 


3(1  +  y*2jl/2  y'y*'2  -  y'**  (1  4.  y’2)2/2  ^  ^  y,2jl/2  yn2 


Simplifying,  one  finds: 


3y'y"2  _  y*»*  ^  4.  yt2j  m  \/ Q  y"2 


(24) 


The  problem  is  reduced  to  solve  equation  (24),  a  third  order  differential 
equation,  the  solution  of  which  gives  y  ■  y  (x) 

Note  that  the  equation  is  homogeneous  and  that  neither  x,  nor  y  appear  in 
it.  To  solve  it,  let  us  first  consider  y'  ■  dy/dx  as  the  variable  and 
y"  «  d  y/dx  as  the  unknown  function.  Therefore,  let 


y*  ■  t,  independent  variable 
y’'«  p,  unknown  function 

One  has: 


y 


Iff 


d  /d 


dx 


J) 


r-p'p 


The  above  equation  is  therefore  of  the  form: 


F(t,p,p')  ■  0 


Substituting,  one  finds: 

3tp^  -  P'P(1  +  t^)  -  (1/C)p^  ■ 
Dividing  by  p,  there  comes: 

3tp  -  p'  -  p't^  -  p/C*  0 
p'  (1  +  t^)  +  p  (1/C  -  3t)  «  0 


or 


or 


p’  _  3t  -  1/C 


Logp 


j  „ 


0 


(25) 


Log  p  ■  3/2  Log  (1  +  t^)  -  1/C  tan"^  t  +  Cj 

going  back  to  the  original  variables,  one  has: 

Log  y"  *  3/2  Log  (1  +  y'^)  -  1/C  tan”^  y'  (  26  ) 

Equation  (26)  is  a  second-order  differential  equation  which  is  again  homo¬ 
geneous,  and  in  which  y  does  not  appear. 

Consider  now  y  as  the  variable  and  y*  «  dy/ dx  as  the  unknown 

function 


hence: 


-  y  -  p 


Substituting  into  equation  (26),  one  finds  now: 

Log  (p  p')  *  3/2  Log  (1  +  p^)  -  1/C  tan  ^  p 
Log  p  +  Log  p'  »  2/3  Log  (1  +  p^)  -  1/C  tan"^  p 


or: 


Log  ^  *  1  Log  (1 

+  p^)  -  Log  p  -  ^ 

,3/2 

-1 

g-(l/C)tan  ^p 

^  P 

:e,  inverting: 

<*)" — hsrr- 

.  ^(l/C)tan-lp  ^ 

(1+  PT 

_^l/Ctan'*p  jp 

(27) 
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to  solve  the  inlegral,  let  us  make  the  change  of  variable  tan~  p  ■  u 
or:  p  ■  tan  u 


2 

ck>  ■  (1  -f  tan  u)  du 

hence:  y  ■  T - e“'^^(l  +  tan^u)  du 

i  (1  +  tan^u) 


Simplifying,  there  comes: 

y 


/  tan  u 


1/2  ~  • 


but 


- — 2~  •  “ 


1  4  tan  u 

hence:  y-Jl^  cosue 
y  «  Jsin  u  e“/^ 


,u/C 


du 


du 


(28) 


Equation  (28)  can  be  integrated  (Burington,  Equation  312): 
(where  a  ■  1/C,  b  »  1) 

gU/C 

y  ■  - - -  (1/C  sin  u  -  cos  u) 

1  ♦  1/C* 

together  with:  p  ■  tan  u 

Going  back  to  the  original  variables,  one  has: 


(2d) 


.u/C 


(1/C  sinu  -  cosu) 


(30) 


1  ♦  1/C 
tan  u 

Let  us  now  differentiate  the  first  equation  (30)  and  substitute 
into  the  second  one: 
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dy  ■  sin  u  du 


therefore: 


cosu 


X  ■  Je'*'  cos  u  du 


.u/C 


(31) 


Equation  (31)  can  be  integrated  (Burington,  Equation  314): 

gU/C 

finding:  x  ■  f  •»  1/C^'  ^  (v  ^  sinu) 

Therefore,  the  parametric  representation  of  the  desired  curve 
is  as  follows: 


.u/C 


XM- — ^ - R  (  (l/C)cosu  +  sinu)  — - — » 

1  +  1/C^  1  + 

-u/C 

y  ■  (  (1/C)sinu  -  cosu)  +  — — «- 

1  +  1/&  1  + 

^-Unu 


(32) 


This  equation  is  plotted  in  Figure  2  for  C>b0.5,  C«1.0,  C«2  and  C  >■  4. 
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m  PROCEDURE  FOR  SOLUTION 


The  equation  to  be  solved  is: 

(A  7/  +  1)  F'2  +  F  F"  +  ^  F'”  ■  0  (14) 

An  approximate  solution  of  this  equation  was  given  in  Reference  5. 
Here,  a  numerical  but  exact  solution  is  discussed,  which  makes  use  of  the 
IBM  7090  digital  computer. 

The  equation  is  solved  using  the  Runge-Kutta  method.  Putting: 

F'  ■  G  (  33  ) 

F”-  G’  -  H  (  34  ) 


equation  (14)  becomes: 

(A  TJ  ♦  1)  G^  +  F  H  -I-  ^  H’  «  0 

or:  H'*  -  2[(A7J  +  1)G^+  FI^  (35) 

The  system  of  equations  (33)  to(3S)  can  then  be  solved  numerically  using 
the  digital  computer,  after  the  initial  boundary  conditions  are  specified. 
A  thorough  discussion  of  the  boundary  conditions  will  now  be  made. 

The  thee  boundary  conditions  suggested  and  used  in  Reference  5 
are  as  follows: 


at 

7J  -0 

F(O).0 

(36) 

at 

TJ.p 

r(o)>G(0)-  1 

(87) 

at 

F’H-GH-O 

(88) 

Preliminary  runs  on  the  computer  indicated  that  the  above  boun* 
dary  conditions  are  neither  correct,  nor  sufficient  to  solve  the  problem. 
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For  these  early  rum,  the  following  procedure  was  adopted:  the  numerical 
calculations  start  from  7^  ■  0.  One  needs  therefore  to  know  the  initial 
conditions  on  the  function  F  and  its  derivatives  at  7^  ■  0,  i.  e. : 

F(0),  G(0),  H(0) 

Equations  (16)  and  (17)  give: 

F(0)  •  0 
G(0)  -  1 

There  ren^ains  to  assume:  H(0)  •  Ct 

From  equation  (7),  if  one  plots  u,  x-component  of  the  velocity,  against)^  , 
H(0)  ■  0(  represents  the  slope  at  7|  «  0  of  that  velocity  distribution. 

For  the  straight  jet,  H(0)  >  0  (the  velocity  profile  is  symmetrical  with 
respect  to  y  ■  0).  For  the  curved  jet,  H(0)  can  have  any  value  between 
0  and  .  For  small  curvatures,  however,  it  stands  to  reason  that 
H  will  be  small.  For  example,  for  A  »  -  0. 25,  one  found  H  zr  -  0. 13. 

One  can  then  start  from  an  arbitrary  value  H,(0)  «  CX,  and 
integrate  the  equations  (33)  to  (35)  numerically.  It  was  found  that  incre¬ 
ments  of  7|  equal  to  0.02  were  satisfactory.  In  principle,  one  inte¬ 
grates  until  71  •  ^  .In  practice,  for  computational  purposes, 

"infinity"  can  be  defined  as  the  value  of  7|  for  which  G(  7|  )  »  0  t  0.000  005. 
It  was  further  found  that,  instead  of  using  a  variable  (  ^  )  for 
which  the  above  criterion  would  be  satisfied,  it  was  just  as  accurate  and 
simpler  to  carry  out  all  calculations  up  to  7|  *  4  and  take  infinity  at 
7|  ■  4.  This  viewpoint  was  carried  out  in  all  subsequent  calculations. 

An  additional  difficulty  which  occurred  because  of  the  use  of  7|  «  4  in 
lieu  of  "infinity"  is  discussed  later  on.  The  trial  values  OC , 
of  0(  can  be  adjusted  until  the  boundary  condition  (38)  at  infinity  is 
satisfied  with  the  desired  accuracy. 
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The  printout  of  early  programs  consisted  of  four  columns  of 
numbers  showing  respectively  7| ,  F,  G  and  H  for  all  values  of  7^ 
between  0  and  4  with  (^7^ «  0.02).  The  velocity  distribution  across  the 
jet  in  each  case  could  be  obtained  by  plotting  G  against  7^  ,  as  shown  in 


It  immediately  became  apparent  that  the  three  boundary  conditions 
(36)  to  (38)  were  not  sufficient,  but  that  one  also  needed  to  have  the 
boundary  condition  H(  »« )  •  0  (i.  e.  the  slope  of  the  velocity  profile  at 
infinity  had  to  be  zero),  in  order  to  have  physically-possible  velocity 
profiles.  Thus,  it  would  appear  at  first  sight  that  one  had  now  four 
boundary  conditions:  F(0)  »  0,  F*(0)  »  1,  G(op  )  «  0,  H(oo  )  ■  0  and  that 
the  problem  would  be  overdetermined.  Such  is  not  really  the  case. 

The  difficulty  is  due  to  the  fact  that  one  had  to  replace  infinity  by  a  finite 
value  of  7|  (7^  «  4),  in  order  to  save  lengthy  step-by-step  integrations. 
It  turns  out  that  the  boundary  condition  G(  )  »  0  is,  for  numerical 
puiposes^equivalent  to  the  two  conditions:  G(4)  «  0,  H(4)  «  0.  If  one 
were  to  be  at  infinity,  G(  oo  )  **  0  would  automatically  entrain  H(  <<>  )  **  0, 
which  therefore  would  not  have  to  be  prescribed  as  part  of  the  problem. 

At  7|  *4,  if  one  were  to  specify  only  one  of  the  two  conditions: 
G(4)  ■  0,  H(4)  ■  0,  either  of  the  two  situations  pictured  hereunder  could 
take  place,  both  being  physically  unacceptable: 
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In  all  subsequent  calculations,  it  was  therefore  accepted  that 
the  boundary  conditions  of  the  problem  at  infinity  could  be  replaced  by 
G(4)  ■  0,  H(4)  «  0.  In  fact,  one  sometimes  had  G(4)  »  0,  but  one  never 
could  have  H(4)  ■  0.  One  usually  had  to  be  satisfied  with  H(4)  4  0.003. 

The  next  step  in  the  work  was  to  set  up  an  iteration  procedure 
for  H(0)  aCX^.  This  could  be  done  simply  by  assuming  two  initial 
values C(, and C<  2  knd  setting  up  an  interpolation  formula  from  which 
the  desired  value  of  CX  could  be  determined.  This  formula  is  not  discussed 
here,  but  was  included  in  the  first  program,  which  is  reproduced  in  part 
IV  of  this  report.  Convergence  was  very  rapid  (usually  2  or  3  iterations). 

During  the  early  runs,  an  alternate  method  of  approach  was 
tried.  The  iteration  was  carried  out  backwards  from  Tf  >4  10  >0. 

The  initial  boundary  conditions  were  then: 

>  assumed  value 

G(4)  -0  ,4 

H(4)  -  0 

The  value  of  was  changed  until  a  value  was  found  which  satisfied  the 
boundary  condition  G(0)  >  1. 

This  method  turned  out  to  be  unsatisfactory  for  many  reasons. 

At  best,  it  had  no  marked  advantage  over  the  direct  integration  method. 
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At  worse,  it  was  very  inaccurate,  because  the  result  d^ended  greatiy 
upon  the  assumed  values  of  G(4)  and  of  H(4).  If  one  assumed  G(4)  ■  H(4)«0, 
the  numerical  integration  would  never  get  "started",  because  these  con¬ 
ditions  could  only  be  satisfied  at  infinity.  One  had  to  assume,  for  example, 
H(4)  =  0.000  005.  Then,  it  turned  out  that  the  whole  solution  depended 
greatly  upon  the  assumed  value  of  H(4).  The  backwards  integration  was 
therefore  abandoned. 

To  summarize  these  early  investigations,  for  any  given  value  of 
a  parameter  A,  the  equation: 

(A  7^  +  1)  F’2  +  F  F"  +  ^  F"’  ■  0 
subject  to  the  boundary  conditions: 

F(0)  -  0 

G(0) «  1 

H(oo  )  •  0 

gives  a  unique  solution:  F(7|),  G(  7^ ),  H(7^  ).  This  solution  is  obtained 
by  determining  by  an  iteration  formula  the  quantity  Qf  ■  H(0),  such  as  to 
satisfy  simultaneously:  G(4)  s::H(4) 

After  a  number  of  runs  were  made,  it  became  apparent  that  the 
above  formulation  of  the  problem  is  not  correct.  Naturally,  in  retrospect, 
this  became  obvious:  it  is  interesting  to  see  how  a  numerical  investigation 
can  lead  to  a  better  theoretical  understanding  of  a  problem!  None  of  the 
statements  made  above  are  wrong;  however,  an  important  boundary 
condition  has  been  omitted;  that  which  expresses  the  matching  on  the 
axis  of  the  "inside"  and  of  the  "outside"  Jet.  For  a  given  value  of  the 
constant  A,  say  A^,  the  inside  of  the  jet  is  found  as  solution  of  the 
equation: 

(Aj  7(1  1)  F'2  ♦  F  F"  ♦  ^  F"  -  0  ; 
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the  outside  of  the  jet  is  found  as  solution  of  the  equation: 

(-Aj  7J  +  1)  F’2  +  F  F”  +  2  ■  0 

These  two  equations  are  obviously  distinct.  If  one  solves  them  together 
with  the  boundary  conditions  (36)  to  (38),  one  must  now  express  in  addition 
that  the  slope  H^(0)  of  the  velocity  distribution  for  the  inside  jet  is  equal 
in  absolute  value  to  the  slope  Hq(0)  of  the  velocity  distribution  for  the 
outside  jet.  Numerical  results  as  well  as  intuition  tell  us  that  there  is 
no  reason  why  this  condition  should  turn  out  to  be  satisfied  automatically; 
the  problem  as  stated  above  is  therefore  overdetermined. 

A  reexamination  of  the  hypotheses  indicated  that  the  problem  can 
be  stated  correctly  for  the  complete  curved  jet  by  eliminating  boundary 
condition  (36):  F(0)«0.  This  condition  was  written  initially  by  analogy 
with  the  straight  jet  case,  in  which  case  it  was  a  consequence  of  symmetry, 
since,  for  the  straight  jet,  v  ■  0  on  the  axis  (there  is  no  cross  flow  across 
the  axis) .  Because  of  the  formal  resemblance  between  the  equations  for 
straight  and  for  curved  jet,  here  again  F(0)  c  0  is  equivalent  to  have  v»  0 
at  »  0,  i.  e.  on  the  axis  (this  can  be  seen  immediately  from  equation 
(8)  ).  This  condition  would  be  justified  if  the  axis  was  a  streamline. 
However,  as  eiqilained  before,  the  axis  of  the  jet  (line  7|  -  0)  is  defined 
by  the  condition  2  K  y  ■  A  7^  ,  and  this  does  not  make  it  a  streamline. 

The  boundary  condition  F(0)  »  0  can  therefore  be  dropped  without  inconsist¬ 
ency  and  replaced  by  the  condition  that 

W®>1  inside  •  ••<“»]  outside  < ’»  > 

<9<ljc0  ®»*jaO 

In  the  remaining  part  of  the  investigation,  the  boundary  conditions 
(37)  to  (3d)  were  therefore  used  in  connection  with  the  numerical  solution 
of  equation  (14).  The  complete  solution  was  obtained  in  two  steps,  which 
will  now  be  described  in  detail. 
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To  obtain  a  numerical  solution  of  the  problem,  It  is  still  necessary 
to  start  with  numerical  values  of  F(0),  G(0)  and  H(0).  As  before,  G(0)  «  1 
and  H(0)  can  be  obtained  by  iteration,  assuming  initially  two  arbitrary  values 
CX  ^  and  Of  2>  would  have  been  possible  to  set  up  an  iteration  proce¬ 
dure  for  F(0)  similar  to  that  developed  for  H(0).  Instead,  it  was  found 
simpler  to  find  the  proper  value  of  F(0)  graphically.  A  first  program  was 
written,  simply  to  establish  in  all  cases  the  proper  initial  values  F(0), 

G(0),  H(0).  The  graphical  procedure  is  illustrated  in  figure  3,  corresponding 
to  the  following  initial  conditions: 

► 

F(0)«  0.001,  0.005,  0.01  successively 
.  G(0).l 

F'(  oo  )  ■  0 

A«+0.25;  A«-0.25 

In  each  of  the  6  cases  above,  equation  (14)  was  integrated  numerically  and 
H(0)  found  by  iteration.  The  result  is  plotted  in  figure  3  in  the  form  of  F(0) 
against  H(0).  In  this,  as  in  all  subsequent  cases,  it  was  found  that  F(0) 
varied  linearly  with  H(0).  The  value  of  F(0)  corresponding  to  the  solution 
of  the  complete  jet  is  found,  according  to  boundary  condition  (39),  as  the 
intersection  of  the  two  lines  A  *  -f  0. 25  and  A  <■  -  0. 25  on  figure  3. 

The  above  calculations  were  repeated  for  a  number  of  values  of  A 
between  0.025  and  0.50.  Some  of  the  results  are  shown  in  figure  4.  Next, 
the  solutions  of  the  curved  jet  problem  F(0)  and  H(0)  found  as  described 
above  were  plotted  against  the  parameter  A  The  results  are  summarized 
in  figure  5.  One  startling  result  came  out  unequivocally:  H(0)  varies 
linearly  with  A;  i.  e.  the  slope  of  the  velocity  distribution  curve  at  ■  0 
varies  linearly  with  the  parameter  A  which  determines  the  dupe  of  the  jet. 

The  variations  of  F(0),  G(0),  H(0)  as  a  function  of  A  are  interesting 
to  follow  on  figure  5,  starting  from  A  ■  0.  A  ■  0  corresponds  to  the  straight 
Jet;  in  this  case,  as  we  know,  F(0)  ■  0  ,  G(0)  ■  1,  H(0)  *0.  As  A  increases. 
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F(0)  becomes  negative  while  decreasing,  then  starts  increasing,  becomes 
zero  and  next  becomes  positive.  This  means  that  the  velocity  v  across 
the  axis  of  the  jet  is  positive  for  small  curvatures  and  becomes  negative 
fur  large  curvatures  (this  comes  out  of  equation  (8)  ).  As  A  increases, 

G(0)  remains  always  equal  to  1  and,  as  noted  before,  H(0)  increases 
linearly  with  A. 

Figure  5  therefore  completely  depicts  all  initial  conditions  of 
the  problem  which  automatically  satisfy  the  boundary  conditions  (37)  to 
(39).  The  information  of  figure  5  can  therefore  be  used  in  lieu  of  the 
boundary  conditions,  hence  reducing  the  machine  program  to  a  straight¬ 
forward  numerical  integration  without  iterations  or  cross-plots. 

A  second  program  was  therefore  prepared  to  perform  the  Runge- 
Kutta  integration,  for  known  initial  values  F(0),  G(0),  H(0).  This  program 
is  listed  in  part  IV  of  this  report.  Note  that  the  information  contained  in 
the  first  program  is  also  listed  in  the  second  one,  but  is  not  used. 

This  was  done  so  that,  if  somebody  wished  at  a  later  date  to  obtain  F(0) 
and  H(0)  with  great  precision,  this  could  be  done  without  having  to  revert 
to  the  first  program.  Therefore,  the  first  program  is  listed  in  part  IV 
for  information  only,  since  all  that  it  contains  is  also  contained  in  the  second 
program. 

Prov^ons  were  made  in  the  second  program  to  calculate  the 
integral  T G^(  7^  )  d7^  and  print  the  result.  This  allows  to  calculate 
the  pressu^  distribution  across  the  jet  as  follows: 

From  equation  13,  an  elemental  pressure  variation  across  the  jet 
is  given  by: 

I  4- 

Remembering  that:  1 "  > 

one  has:  dp  ■  ^  PU^  A  G^  d 7^ 
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The  total  pressure  across  the  jet  is  therefore: 


Ap-j  P  u*  a 


Gh7l)dri 


(40) 
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IV  PRCX3RAM 


Two  programs  are  showa  here.  The  first  program  was  used  to 
obtain  the  data  of  figure  5.  The  second  program  is  more  general  and 
includes  the  elements  of  the  first  program  as  an  option. 

All  symbols  used  are  defined;  listings  of  the  two  programs  are 
shown,  as  well  as  the  data  card  input. 


Symbolic  location 


Contents 


A 

ABSF 

AF 

ALP  ( 7;  ) 
DEL  X 
EN 

Fill) 

F  1 

G(T|) 

G  1 
H(7|  ) 

II 

IIJ 

HP 

INT 

INTM 

lOP 

JJ 

K 

M 


Constant  in  differential  equation  (12) 

Closed  subroutine  used  to  compute  the  absolute  value 
Initial  value  of  F  (i.  e.  F(0)  ) 

Initial  value  of  H  (i.  e.  H(0) ) 

Increment  of  7|  (DEL  X  ■  0.02) 

Maximum  numerical  value  of  7^  (EN  =  4) 

Unknown  function  of  the  problem 
.Tl  at  specified  interval  of  integration 
-  F’  (/?  ) 

Value  of  G  from  the  previous  integration  step 

{71) 

index  of  values  of  A 
iteration  counter 

index  variable  for  the  value  of  F(0) 

closed  subroutine  which  initiates  the  subroutine  RWINT 
(SHARE  distribution) 

subsequent  entries  to  RWINT 

Note:  DAUX  is  used  by  RWINT  to  compute  the  deriv¬ 
atives  G  and  H. 

print  option:  if  non  zero  print  every  case 

if  zero,  only  print  the  converged  value  of  H 

variable  index  used  to  go  from  one  value  of  H  to  the  next 
one 

index  variable  used  in  interpolation  formula 
the  number  values  of  A 
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Symbolic  location 


Contents 


MM 

Maximum  number  of  iterations  and  maximum  number 
of  values  read  in  for  F(0) 

OP 

option 

T(2) 

n 

T(3) 

delta  7^ 

T(4) 

F(0) 

T(5) 

G(0) 

T(6) 

H(P) 

Additional  symbols  for  the  second  program 


G  11 
G  12 
NN 

SIMPUN 


G  square 

sum  from  ^,to  7^ g  ^  square  d 7^ 
option  whether  to  compute  G  12  or  not 

2 

Closed  subroutine  to  evaluate  the  integral  sum  of  G  d  /7 
(SHARE  distribution). 
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HRST  PROGRAM 


DIMEMSION  A(  JOO)  ,T(  100)  ,Ft  1000)  ,i;(  1000)  ,H(  1  000  )  ,  ALP  1 1  1  00  )  ,  ALP  (  1  000 
1)  ,Fl(b00) 

X.AFdOO) 

COHMOi^  A,T,F,li,H,I  ItALPl.ALO 
X,F  1 

16  FORMAT  (78H0  FTA  F  G 

X  H  ) 

READ  10 

10  FORMAT  {72H1 

X  ) 

PRINT  10 

READ  1?,EN,M,MM,DELX, ( ALP( JJ) , JJ=1 f 2) , lOP 
12  FORMAT  ( 1F6.C,2I6,3F16.6,I6) 

READ  15,  (All ) , 1=1 ,M) 

15  FORMAT  (UF16.6) 

READ  1 70, ( AF( I ) , 1= 1 ,MM) 

170  F0RMAr(4El6.4) 

PRINT  50,£M,M,MM, ( ALP( JJ) , JJ= 1 ,2) ,DELX 

50  FORMAT  (25H0  THE  UPPER  LIMIT  ON  ETA=F6.2/27H0  THE  NUMBER  OF  VALUE 
XOF  A=I6/50H0MAXIMUM  NUMBER  OF  I TERAT  I0NS= I6/22H0  UITIAL  VALUES  OF 
X  H=2F 16.6/)2H0  STEP-S I ZE  =  F 1 6. 6 ) 

PRINT  51 

51  FORMAT  ( 1H0,50X,8HA-VALUES  ) 

PRINT  52, (Ad), 1=1, M) 

52  FORMAT  (50X,F16.6) 

I  1=1 

IFCIOP)  53,55,55 
55  SENSE  LIGHT  2 

53  SENSE  LIGHT  1 
J  J=  1 

I  IP=1 
K=1 
IIJ  =  1 

IF  (EN)  14,14,15 

14  CALL  ENDJOB 

15  T(2)=0.0 
T (5)=DELX 
T(4)=AF( IIP) 

T(5)=1.0 
T(6)=ALP( JJ) 

1  =  1 

CALL  INT(T,3, 1 ,0,0,0,0,0,0) 

F1(I)=T(2) 

F(I)=T(4) 

G( I )=T(5) 

H( I )=T(6) 

1  =  1  +  1 

77  CALL  INTM 
F1(I)=T(2) 

FII)=T(4) 

G( I )=T<5) 

19  H(I)=T(6) 

1*1  +  1 

IF(EN-T(2) ) 18,77,77 
18  IF(ABSF(T(5))-.0005)  59,39,40 

39  SENSE  LIGHT  0 
GO  TO  22 

40  IFtSENSE  LIGHT  1)  41,45 

41  JJ*JJ+1 
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GI=T(5) 

IF  (SENSfc  LIGHT  2)  21, IS 

45  ALP(  K*2)*  {ALP(K*n-(T(5)»(ALP(K)-ALP(K  +  1  )  )/(Gl-T(  5))  )  ) 
JJ=K+2 
K*K+I 
G1»T(5) 

I IJ-IIJ+1 

IF(IIJ-MM)  56,56,80 
56  IFISENSE  LIGHT  2)  21,15 
80  PRINT  25,IIJ 

25  FORMAT  (20X,I6,21H  ITERATIONS  COMPLETED) 

SENSE  LIGHT  0 
I  IJ=1 
GO  TO  22 

21  SENSE  LIGHT  2 

22  PRINT  10 
PRINT  16 
PRINT  17,A(II) 

17  FORMAT  (  1HO,40X,F16.6) 

1  =  1 

23  PRINT  31,F1(I),F(I),G(n,H(I) 

IFCEN-Fl (I) 130,30,32 

30  IFISENSE  LIGHT  2)33,35 

32  I  =  U) 

GO  TO  23 

33  CONTINUE 
SENSE  LIGHT  2 
IF(IIJ-MM)  15,15,35 

35  IIP*1IP*1 
I  IJ=1 

IF(IIP-MM)  190,180,180 

190  SENSE  LIGHT  1 
SENSE  LIGHT  2 
K») 

JJ=1 

GO  TO  15 
180  IIP»1 
I  IJ  =  1 
11=11+1 

IF(A(II))  191,192,192 

191  JJ=1 
ALP(JJ)=-ALP{JJ) 

ALP(  JJ  +  1  )=-ALT»(  JJ+1 ) 

GO  TO  193 

192  JJ=1 

ALP( JJ)»ABSF{ALP( JJ) ) 

ALP ( JJ+1 1»ABSF( ALP (JJ+1 ) ) 

193  SENSE  LIGHT  1 
SENSE  LIGHT  2 
K-1 

JJ*1 

IF(II-M)15,I5,14 

31  FORMAT  (4F16.6) 

END 
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SLOROUTINE  UAOX 

DI><CNSICN  A(  IOO)tT(  ICO)  t>^(  1000)  i'J(  1000)  fH(  1000)  lALPI  (  100)  .ALP(  1000 
1)  ,H(500) 

COMMON  A,T,F,G,H,I1,ALPI,ALP 
TI7)=Tt5l 
1 (d)=r(6) 

T(Vl»-2.( (A(II)*T(2)*1. )*TC5)»T(5)+T(4I«T(6)) 

KETUKN 

ENu 
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SECOND  PROGRAM 


f.'IfibfMSIO  .  A(  100  )  ,T  (  100)  ,F<  1000)  i  ,(  I  0  00  )  ,  H(  100  0  )  ,  ALP  I  (  1  00  )  ,  (  1  000 

1)  .DCbOO) 

X,AP( 100) 

X,  01  1  (  bOC)  ,G12rjC0) 

CLMMlM  A  ,  f  ,  ,  0,H,  I  I  .AL-^  I  ,  AlP 
X,  F  I 

X,ol  l,GI? 

16  FCaPAT ( IPOHO  FTA  F  G 

X  ('  GtSfJuA,<Ei  SUM!  G(  (S  )UA^E)  XO  FFA  )  ) 

RtAF  10 

10  FORMAT  (72H1 

X  ) 

PRINT  10 

REAP  1?1  K'P,  J'J.OLLX 

121  F0RMAT(i*16,F16.6) 

REAL  12,  EIM,(ALPtJJ),JJ=),IOP) 

12  F(.RFAT  ('IF16.6) 

REAL  H,  (All), 1=1, M) 

13  FORMAT  (  t<F  16.6) 

READ  170, {AF(  I  ) ,  '=  1,HM) 

170  FORMAT  (4tl6.l») 

PRINT  60,EN,M,MM,(ALP(JJ),JJ=1,2),0£LX 
bO  FORMAT  ( 25H0  THE  UPPER  LIMIT  ON  c:TA  =  Fo.  2/2  7HD  The  \iuM3cR  OF  VALUE 
XUF  A=I6/30HflMAX  IMUM  NOF'RtR  OF  1 1  £R  A  T  1 0\ S  =  1  6  /  2  ?H  0  INITIAL  VALUES  OF 
X  H=2F16. 6/121.0  STEP-SI/E^FlV-.o) 

PRINT  51 

51  FORMAT  ( IHC,30X,6HA-VALUES  ) 

PRINT  52,  (  A(  I  )  ,  1  =  1  ,M:) 

52  FORMAT  ( .•>0X,F  16.6) 

11=1 

IFIIC'P)  •j3,SA,b5 
55  sense  light  2 

53  SINSE  LIGHl  1 
JJ=  1 

I  IP=1 
K=  1 
1  I  J=1 

IF  (EN)  14,  111,  15 

14  CALL  EnOjOH 

15  T (2  1=0.0 
Tt3)=OELX 

T {4)=AF( IIP) 

T ( 5)= 1 .0 
T (6)  =  ALP( JJ  1 
1  =  1 

CALL  INTIT, 3, 1,0,0,0,0,0,01 

r  )  I  I )  =  r  (  2- ) 

F I  I )  =  T(U) 

G( 1)=TI5) 

HI  I )  =  TI6) 

G I  1  I  I )=3I I )*CI I  ) 

1  =  1-*  ) 

77  CALL  INTM 
F1I I)=T|2) 

FII)=TIU) 

Gl 1 )=TI5) 

19  h: I)=TI6) 

Gil  I  1 )=GI I  )*&l I ) 

1  =  1*1 
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IF( tN-TI2 ) ) ItOO, /7,  {7 
IbOO  IF(N'J)  Ih, lb, 1801 
IbOl  NrMl=2 

Du  1832.  l  =  l,NfJ 

IH02  GI2(  l)=SiyPUiM(H(  I  )  ,G1  1  (  n  f  N\1  ) 

PiliNr  10 
PIMNT  16 
00  180^  1=1, GU 

idOi  PKii^r  1  il  ,n  ( I  ),F{  I ) , ;( I )  ,H(  I )  ,G  1 1(  I )  ,r,i?(  I ) 

Id  I  FCUMAr(6F 16.61 
JJ  =  JJ+  1 

I  1.^=1  IP*  1 

II  =  IU1 

IF( JJ-IOP) IS, IS, 14 
18  IF(fi8SF(  I  (  j)  )-.000a)  59,39,143 
.39  St.MSfe  LIGHT  0 
GO  TO  22 

40  IFISE'JSE  LIGHT  1)  4l,4S 

41  JJ=JJ+I 
G  1  =  r  (  S  1 

IF  (ScHSE  LIGHT  2)  21,  IS 

4S  AlP(  K  +  2)=  (/'LPtKH  1 )-(  T  (S)  •(  ALP(K)-ALi'(  <+l  )  ) /(  Gl-T(  5)  )  )  ) 
J J=K+2 
K  =  K+1 
G1=T(S) 

I 1J=I IJ+1 

IF(I1J-I‘y)  S6,S6,8G 
56  IFlSi.'ISE  LIGHT  2)  21,  IS 
SIC  PRINT  25,  IIJ 

25  FORt-.AT  {?3X,I6,2IH  11E:UTIUNS  CONPLETllD) 

SENSE  LIGHT  0 
I  IJ=  1 
Go  TO  22 

21  ScMSE  light  2 

22  PRl a  10 
PRI.mT  16 
PRI^iT  17,  A(  11  ) 

17  FORNAT  ( 1HC,40X,F 10.6) 

1  =  1 

23  PRINT  31 ,F 1 ( I ) ,F ( I ) ,G( I ) ,H( I ) 

I F( E\-F1 ( I ) )3T, 3C, 3^ 

30  IFISENSE  LiGi-T  2)33,35 

32  1=1+1 

GO  TO  23 

33  CONTI XJE 
sense  Lli.HT  2 
IF(IIJ-MI-.)  IS,  15,3s 

35  I IP=I IP+  I 
I  IJ=  I 

IF(IIP-MT.)  190,1uC,l80 
190  SENSE  LIGHT  I 

IF(IOP)  1900,1903,190 
1901  SLNSL  LH.HT  2 
1900  K=1 
JJ=  1 

GO  TO  15 
180  IIP=I 
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1  IJ=  I 
11=11+1 

!F(A(II))  I91,19/!,1V2 

191  JJ=1 
ALF(JJI=-ALP( JJ) 

ALH(.1J  +  1  )=-ALP(  J  )+  1  ) 

no  10  19^ 

192  JJ=I 

ALP(JJ)  =  AHSF  (ALP(JJ)  ) 

ALP(  JJ  +  n=APSF(  ALPt  JJ  +  I  )  ) 
19,1  SbNSE  LIOHf  1 

IF(IOP)  1914,1931,19^0 
1930  S£i\S£  LIGHT  2 
19.11  K=l 
JJ=  I 

IF(  1  I-M)  IS,  !•),  14 
11  FOKMAT  (4F16.6) 

END 
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Input 


For  the  second  program,  the  input  data  consists  of  the  followii^ 

cards: 


Card  No. 

Columns 

Contents 

1 

1 

1 

2-72 

Problem  identification 

2 

1-6 

M 

7-12 

MM 

13-18 

lOP 

18-24 

NN 

25-40 

DELX 

3  to  N 

1-16 

EN 

i, 

1 

17-NN* 

H{0) 

1  to  M 

1-NN 

values  of  A 

M<f  1  to  L 

1-NN 

values  of  AF 

*  The  first  value  of  H(0)  is  punched  in  Columns  17-32,  and  successive 
lOP  values  are  punched;  field  width  of  16  up  to  and  including  column  NN; 
4  words  per  card. 


For  Che  first  program: 

same  as  the  second  program,  except 

lOP-  2 
NN  *0 

leave  cards  N  -f  1  to  L  off. 
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V  DISCUSSION  OF  RESULTS 


As  discussed  earlier,  the  emphasis  in  the  present  work  was 
placed  more  on  the  formulation  and  the  solution  of  the  problem  than  on 
the  exploitation  of  the  results.  The  program  included  in  part  IV  of  the 
report  will  allow  anyone  to  generate  immediately  as  many  velocity  and 
pressure  distributions  as  he  wishes. 

Typical  velocity  and  pressure  distributions  are  shown  in 
Figures  6  and  7  for  several  values  of  the  curvature  parameter  A.  The 
first  figure  shows  a  plot  of  G(  7^  )  against  7^  .  This  only  represents 
the  x-component  of  the  velocity;  in  order  to  have  the  true  velocity 
distribution,  it  would  be  necessary  to  add  the  y-component  of  the 
velocity.  The  plot  is  representative,  however,  since  v  is  eiqpected 
to  be  small  compared  to  u.  The  velocity  distribution  plot  shows  that, 
as  the  curvature  increases,  the  maximum  velocity  U  increases  and  the 
location  of  the  point  of  maximum  velocity  moves  radially  outwards  (in 
the  region  of  positive  curvature,  or,  for  an  annular  jet,  towards  the 
center  of  the  cavity).  This  trend  is  in  agreement  with  the  experiment¬ 
al  data  on  two-dimensional  curved  jets  of  Reference  7.  Overall, 
however,  the  differences  between  straight  and  curved  jet  velocity 
distributions  are  not  substantial  (again,  this  could  be  modified  due  to 
the  effect  of  the  transverse  velocity). 

Figure  7  shows  a  plot  of  the  calculated  pressure  distribution 
across  the  jet.  It  can  be  seen  that  the  pressure  variation  is  not 
exactly  linear  as  the  momentum  theory  would  predict.  Also,  the 
pressure  rise  on  both  sides  of  the  axis  of  the  jet  are  not  exactly 
equal. 
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laCBeaBBBBfwBaUBBBaBBBBaBBBBBBBBBaBBBBBBBBBBBaBBBBBaflBBBBBBBaaBBBBBB' 
JBBBBBBBBBWaBaBaBBBBEaaBBBiBBBaBaBaaBBBBBBBBBBBBBiBBBBBaBBBBBBBBiBBB 
laBBIBBBBBBBiriiaBBBBaaBaBBBBBBBBBBBBBBBBBBBaBBBBBBBBBBBBBiBBBBBBBiBBBBB 

|aB*'B.LaBBBiaa».'iaBaBBBBBaBBBaBBaBBBBBBBBaJaBBBBBBBBaBBBBBBBBBBBBBBBBaBBB 
iBBiWfilBBBBBBJIBBBBBBBBBBBBBBBBBBaBBBBBBa.iaBaBBBBBBBBBBBBaBBBBBrBBaaBBBB 
laBMB^BBaBaaB■aaaBaBaaBaBaaBaBBBBaBBBaBB^BBBBBBaBBBBaBBBBBBBBB^- BBBBBBBB 
IbB^B I BBBBBBBBBBBBBBBaBBaBaBBaaBaBBBaBBBV  BBBBBBBBBBBBaBBBBBflB  BalBBBBaBB 
laB'<BaBBBBBBBi;aBaBBaaaBBBBBBBaBaBaBBBBBBKBBBBBBBBBBBBBBBBBBBBBi:jBBBBBBB 
iBBiiBBBBBBBBBBBBBBaBBBBBBBBBBBBBBBBiBBBBnBBBBBBBBaBBaBBaBBBBBBaBBBBBBBa 
larr.aiBBBaBBa.jBaaBaBBBto="BBBBBBBBBaBBBBBrBBBaBBBBBBBaBBBBBBBB  BBBBBBBBBB 

SB  mm iBBBBBBr  ;BBBBBBaaBBBBk=?BaaBBBBBBBV~«BBBBBBBBBB9BBBBBBBB  BBBBBBBBBB 
■  '.a  r  IBBBBaBkBBaBaBBaBBBBBBBaB»C?BBaBBaL*IIBaBBBBBBBBBBBBBBBBB  BBBBBBBBBB 
|Brf«-.BaBaBBBBBBBBaBaBBBBBBBBBaPBBBB^=^BBB?BBBBBBBBBBBBaBBBBBBiBBBBBBBB 
iBrUBSnBBBBBBBBBaBBBBaBaaBaBBBB IBBBBBBBBBf^^BBBBBBBBBBBBBBBBBBBBBBBBBBB 
iBaBaJBBBBBBBBBBBaBBBBBBBBBBBBBCBBBBBBBBBBBBBkZ^PBflBaBBBBBaBB BBBBBBBBBB 
laBBBBaaBaBaBBBBaBaaBBBBBBaBaBB^  BaBBBBBBaBBBBBBBBMii:=«aBBBBBBB B BBaBBaBB 

ISBaBBBBBBaBBaBBBBaBflBBBBBBBBBBK'BBBBBBBBaBBBBaBBBflBBBBBBiBSh"* 'BBBBBBBB 
BB'i  BBaBBBBBaBBBBBBaBaBBBaBBBBanBBBBBBBBBBBBBBBBBBBaaBBaBBBBBBBB»^-"3BBB 
Bak.BBBBBBBaBBBBaBBBaBBaBBaBBaBiiflflBBBBBBBBBBBBBBBBBBBBBBBBBBP'SSg 
BBBBBBBBBBBBBBBBBBBBBBBaBaBaaanBBBBBBaaBaBBBBBBBBBBaBB»!:=0BBBBi 

aBr'iaBBBBBaaBBBaBBaBaBBBBaaBar’'BBBBBBBBBBBBBBBBaBP;:Sa«aBBBBBBBI. _ 

BB.VBaBBBBBBaBaBBaBBBBBBBBBBBaBL^BBaBBBBBBBBBBP'SaBBBBBBBBBaBB  BBBBBBBBBB 

.BB.aBBBBBBaBBBBBBBaBBBBBBaBBBaBraBaBBBBPn^MBBBaBBaBBBBBBBBBBBBBBBBiBBB 

|BB  I  IBBBBBBBBBBBBBBBBBBBBBBBBBBBByilBP^^aBBBBBBBBBBBBBBBBBBBBBB  BBBBBlEiBB 

_ IBBBBaaaBaaBBBBaaaBBBBMBBT'BBBBBBBBBniBBaBBBaaBnBBaBaBBBaa.rBBBBBBBBl 

|aBBBBBBaBr.>''aaiBaBBBrrnBaBBBBaBi:''aBBBBBBBW  ^BBBBBBBBTlIBBflBBBBBBBBBBBBBBBi 
I  BBBBBBaBBB  :^aB|BBBBai^  BBBB^  JBBBr  ^BBBBBBBBr  .’'BBBBBBBBblsiBBBBBBBBI  il'BaiBBBBBB 
BBBBBBBBBB  •^BBiBBBBBk«'JBBBPaaBBk*.BBBBBBBBL'^BBaBaBaait)IBBBBaaBBk'JBaBBiBBBB 
BaBaBBBBBI'liaaaBBBBardlBaa.V'JBBBl'iaBBBBBBBM'BBBBBBBBr.lBiBBBBBaBaBBBBBBBBB 
BBBBBiBBBB«<BBaBBBBBBraaw-SiaBBB^<aBBBBBBBB<<BBBBBBBaL"iBiBBaBaBBBBBBBBBBBB 
.BBBBBBaBBaBBBBBBaBBBBBBnBBBBaBBBBaBBBBBBBBBBBaBBBB^^BBBBBBBBaBaBBBBBBBBl 
iBBBBBBaBBBaBBaBBBBiBBBBIII  IBBBBBBBBBBBBBBBBBBBBaBBBBaBBBBaBBBBBBBBBBBBBBl 

|aBBaaBiBBaBaBBBaaBiBaBBL;«iiBBBBiBliaBBBBBBaBBBaBBaBBBBBBBBBBBBBBBBBr - 

iBBBBBBiBBBBBBBaaaBBBaaBBBBBBallflEiBBBBBBBBaaBBaBBBBBBBaaBBBBaBaBBBI 
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